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New results for the strength of the symmetry energy are presented which illustrate 
the complementary aspects encountered in reactions probing nuclear densities below 
and above saturation. A systematic study of isotopic effects in spectator fragmenta- 
tion was performed at the ALADIN spectrometer with 124 Sn primary and 107 Sn and 
124 La secondary beams of 600 MeV/nucleon incident energy. The analysis within the 
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Statistical Fragmentation Model shows that the symmetry-term coefficient entering the 
liquid-drop description of the emerging fragments decreases significantly as the multi- 
plicity of fragments and light particles from the disintegration of the produced spectator 
systems increases. Higher densities were probed in the FOPI/LAND study of nucleon 
and light-particle flows in central and mid-peripheral collisions of 197 Au+ 197 Au nuclei 
at 400 MeV / nucleon incident energy. From the comparison of the measured neutron and 
hydrogen squeeze-out ratios with predictions of the UrQMD model a moderately soft 
symmetry term with a density dependence of the potential term proportional to (p/po) 7 
with 7 = 0.9 ± 0.3 is favored. 



1. Introduction 

Microscopic calculations are rather consistent in their predictions for the strength 
of the symmetry energy at nuclear densities below saturation.^ For homogeneous 
matter, the symmetry energy decreases with decreasing density, roughly in pro- 
portion to (p/po) 7 where po represents the nuclear density at saturation and 7 is 
approximately 2/3. The tendency of nuclear matter to cluster modifies this result, 
however. The symmetry energy is ex pected to be finite at very low densities, mainly 
because of a-particle formationpEl anc [ should reflect the abundance of larger frag- 
ments at the freeze-out conditions of multifragmentation reactions. In the Statistical 
Multifragmentation Model which has been successfully applied to many types of 
fragmentation reactions, this is accounted for by using, in its standard version, the 
symmetry-term coefficient of Weizsacker's formula for stable nuclei for the descrip- 
tion of the emerging fragments.^ Recent experiments have shown, however, that 
this scenario of an idealized freeze-out will have to be modified in order to be con- 
sistent with isotope distributions and other observables derived from isotopically 
resolved fragment yields (for references, see reviewi^^ and contributions to this 
proceedings). The excited fragments seem to be modified in the hot environment 
which has important consequences for astrophysical processes such as the collapse 
and explosion of massive stars.^ 

At higher-than-normal densities, on the other hand, ou r kn owledge regarding 
the strength of the symmetry energy is much less complete Most of the avail- 
able results are obtained by extrapolating from below or near saturation to higher 
densities, thereby relying on a power law parameterization (p/po) 7 for the density 
dependence of the potential part of the symmetry term. According to present find- 
ings, the density dependence is moderately soft, with 7 in the r ange of 0.5 to 1.0, 

IQIQI 

and appears to be rather consistent among different observables i 21 ^ Very recently, 
a super-soft symmetry term was deduced from the analysiJ^ of the charged-pion 
ratios for 197 Au+ 197 Au reactions at 400 MeV/nucleon incident energy measured 
by the FOPI collaboration^^ This result is of particu lar int erest because of its 
consequences for neutron-star properties (see, e.g., RefsP^HSI) brtt obviously con- 
tradicts the extrapolations. To clarify the situation, it will be necessary to employ 
additional probes sensitive to the high-density stage of the reaction and capable 
of constraining the model analyses and the deduced parameters of the equation of 
state of neutron-rich matter. 
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The Symmetry Energy in Nuclear Reactions 3 

The new results presented in the following, obtained from experiments con- 
ducted at the GSI laboratory with beams of several hundreds of MeV per nucleon, 
address these currently discussed aspects. From an isoscaling analysis of spectator 
fragmentation, a considerable modification of the symmetry term coefficient rep- 
resenting the fragment properties in the hot freeze-out environment is deduced, 
and the ratio of neutron versus hydrogen squeeze-outs is shown to be a sensitive 
observable for determining the symmetry energy at supra-saturation densities. 

2. Isoscaling in Spectator Fragmentation 




Fig. 1. Inclusive mass distributions for fragments with 3 < Z < 6 obtained for 107 Sn (left panels) 
and 124 Sn projectiles (right panels). 



The ALADIN experiment S254, conducted in 2003 at the SIS heavy-ion syn- 
chrotron, was ai med at studying isotopic effects in projectile fragmentation at rel- 
ativistic energies Q3ESI Besides stable 124 Sn beams, neutron-poor secondary 107 Sn 
and 1 La beams were used in order to explore a wide range of isotopic compositions. 
The radioactive beams were produced at the fragment separator FR£p2 by the frag- 
mentation of primary 142 Nd projectiles with energies near 900 MeV/nucleon in a 
thick beryllium target. The FRS was set to select 12 La and, in the second part of 
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the experiment, 107 Sn projectiles which were then delivered to the experiment. All 
three beams had a laboratory energy of 600 MeV/nucleon and were directed onto 
nat Sn targets with an areal density of 500 mg/cm 2 . The acceptance for fragments 
from the projectile-spectator decays was 90% for fragments with atomic number 
Z = 3 and increased gradually to about 95% for Z > 6. The obtained mass resolu- 
tion was 7% (FWHM) for Z < 3, decreasing to 3% for Z > 6, so that masses for 
Z < 10 were individually resolved (Fig. [1]). 
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Fig. 2. Experimental data (stars) and SMM ensemble calculations (open symbols) of isoscaling 
coefficients a extracted from yield ratios of fragments (3 < Z < 10) from 124 Sn and 107 Sn 
projectiles (top panel), and from 124 Sn and 124 La projectiles (bottom), versus Abound normalized 
with respect to the nominal projectile charges Zq = 50 and 57. Four symmetry-term coefficients 7 
were used in the SMM calculations as indicated in the figure. 

In order to reach the necessary beam intensity of about 1000 particles/s with 
the smallest possible mass-to-charge ratio A/Z, it was found necessary to accept 
a distribution of neighbouring nuclides together with the requested 124 La or 107 Sn 
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isotopes. Their mean compositions were <Z> = 56.8 (49.7) and <A/Z> = 2.19 
(2.16) for the nominal 124 La ( 107 Sn) beams, respectively!^ Model studies consis- 
tently predict that these <A/Z> values are also representative for the spectator 
systems emerging after the initial cascade stage of the reaction. In particular, the 
differences in <A/Z> between the neutron-ric h and neutron-poor cases are ex- 
pected to remain the same within a few percent [HUH] 

Global fragment ation o bservables were found to depend only weakly on the 
isotopic composition! ^^ This includes, in particular, the mean multiplicity of 
intermediate-mass fragments (3 < Z < 20), the largest atomic number Z max within 
a partition, and the evolution of these quantities with Abound- The sorting variable 
■Abound = ^Zi of fragments with Zi > 2 represents the atomic number Z of the 
spectator system, apart from emitted hydrogen isotopes, and is inversely correlated 
with the transferred excitation energy. 
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Fig. 3. Measured cross sections for the fragment production in selected Abound intervals of 10-unit 
width as a function of the fragment Z (symbols) in comparison with the results of SMM ensemble 
calculations (lines) for the case of 124 Sn projectiles. The centers of the five Abound intervals and 
the scaling factors used for displaying the cross sections are indicated in the figure. 



The comparison of the measured fragment yields from neutron-rich and neutron- 
poor systems shows that isoscaling is observed!^ The isoscaling parameter a, de- 
termined from the yields for 3 < Z < 10 is found to decrease rapidly as the 
disintegration of the spectator systems into fragments and light particles increases 
(Fig. [2]), confirming earlier results for the fragmentation of target spectators in re- 
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actions of 12 C on 112 ' 124 Sn at 300 and 600 MeV/nucleonP^ Nearly identical results 
are obtained for the isotopic and isobaric pairs of reactions. 

The analysis of the data with the Statistical Multifragmentation ModeP was 
performed with ensemble calculations adapted to the participant-spectator scenario 
at relativistic energies!^ The ensembles of excited systems with varying excitation 
energy and mass were chosen so as to best reproduce the charge spectra and corre- 
lations observed for the fragment production (Fig. [3]) . The mean neutron-to-proton 
ratios <N>/Z for intermediate-mass fragments up to Z = 10 and the isoscaling 
parameters were found to be particularly sensitive to the coefficient 7 of the sym- 
metry term E sym (A, Z)=MA — 2Z) 2 /A in the liquid-drop description of excited 
fragments at freeze-outP2l23] p or ^ ne j SO scaling parameter a, this is demonstrated 
in Fig. [2j The SMM standard value 7 = 25 MeV is applicable only in the bin 
of largest Abound- Smaller values have to be chosen for reproducing the rapidly 
decreasing parameter a in the fragmentation regime at smaller Abound- 




Fig. 4. Symmetry term coefficient 7 as expected from the changing fraemcnt-mass distributions 
using surface and volume symmetry-term coefficients from Rcfs. ^ ?° (open symbols) in com- 
parison with the coefficient 7 for hot fragments obtained from the isoscaling analysis for the 
107, 124 g n p a j r Q £ reac tions with the Statistical Multifragmentation Model (dots). 



The symmetry-term coefficient 7 for hot fragments resulting from the isoscal- 
ing analysis for the 107,124 Sn pair of reactions is shown in Fig. |4] as a function of 
■^bound/^u- F° r comparison, three predictions are shown in the figure, obtained with 
rather different approaches but all containing the effect of the surface-symmetry 
term whose importance increases for the lower-mass fragments. The coefficients of 
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the mass formula of Myers and Swiatecki are adapted to ground-state masses!^ 
From the energies of isobaric analog states, a relation between the volume and 
surface capacitances of nuclei for absorbing asymmetry N — Z was derived by 
Danielewicz and LeeJ^ while Kolomietz and Sanzhur have used a variational ap- 
proach using Skyrme forces to derive equilibrium values for the volume symmetry 
term with surface and curvature corrections for nuclei along the p-stability lineH^ 1 
With these coefficients, the effective symmetry energy averaged over the set of 
partitions was calculated for the five bins in Abound after the experimental Z dis- 
tributions had been converted to mass distributions using the projectile N/Z. The 
obtained results show similar trends. The smaller fragments produced at higher 
excitations cause the effective mean symmetry term to decrease with decreasing 
inbound in all three cases but at a slower rate than that resulting from the isoscaling 
analysis of the experimental yield ratios. 

3. Neutron and Hydrogen Flows 
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Fig. 5. Elliptic flow parameter V2 for mid-peripheral 197 Au+ 197 Au collisions at 400 MeV per 
nucleon as calculated with the UrQMD model for protons (circles), neutrons (triangles), and 
the total hydrogen yield (stars) as a function of the laboratory rapidity t/i a b- The results have 
been filtered to correspond to the geometrical acceptance of the LAND setup used in the joined 
experiment. The predictions obtained with a stiff and a soft density dependence of the symmetry 
term are given in the upper and lower panels, respectively. 



In two experiments at GSI combining the LAND and FOPI (Phase 1) detectors, 
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neutron and hydrogen collective flow observables f rom 19 7 Au+ 197 Au collisions at 
400, 600, and 800 MeV/nucleon have been measured P^^ This data set is presently 
being reanalyzed in order to determine optimum conditions for a dedicated new 
experiment but also with the aim to produce constraints for the symmetry energy 
at the high densities probed in central collisions at these energies. The results 
reported here are obtained by comparing with predictions of the UrQMD model 
which has recently been adapted to heavy ion reactions at intermediate energies!^ 
The predictions obtained for the elliptic flow of neutrons, protons, and hydrogen 
yields for 197 Au+ 197 Au at 400 A MeV are shown in Fig. [5l Two values are chosen 
for the power-law exponent describing the density dependence of the potential part 
of the symmetry energy, 7 = 1.5 (asy-stiff) and 7 = 0.5 (asy-soft). The UrQMD 
outputs have been filtered in order to correspond to the geometrical acceptance 
of the FOPI/LAND experiment. This produces the asymmetry of V2 with respect 
to mid-rapidity yi a b = 0.448 because higher transverse momenta are selected with 
increasing rapidity. The neutron squeeze-out is significantly larger (larger absolute 
value of v 2 ) in the asy-stiff case (upper panel) than in the asy-soft case (lower 
panel) while the proton and hydrogen flows respond only weakly to the variation 
of 7 within the chosen interval. 




Fig. 6. Directed-flow parameters i>i for neutrons (triangles) and hydrogen isotopes (stars) for mid- 
peripheral (5.5 < b < 7.5 fm) collisions of 197 Au+ 197 Au at 400 MeV/nucleon as a function of the 
scaled rapidity y/y p in comparison with the UrQMD predictions for 7 = 1.5 (a-stiff) and 7 = 0.5 
(a-soft) represented by the dashed lines. 



A nearly negligible sensitivity to the stiffness of the symmetry energy is ex- 
hibited by the directed flow, according to the UrQMD model. In Fig. [fU the first 
Fourier coefficient v\ of the azimuthal particle distribution with respect to the re- 
constructed reaction plane is shown as a function of the rapidity y, normalized 
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with respect to the projectile rapidity yo = 0.896, and for the interval of impact 
parameters 5.5 fm < b < 7.5 fm. The dashed lines representing the results for 
neutrons and for hydrogens and for asy-soft and asy-stiff parameterizations of the 
symmetry energy fall practically on top of each other. The predictions compare well , 
however, with the experimental results for the multiplicity bin PM3 (see RefP3) 
expected to correspond to this range of mid-peripheral impact parameters (Fig. [6]). 
The range of rapidities y > 0.8 at which the data deviate from the expected lin- 
earity coincides with the shift of the LAND acceptance to transverse momenta 
Pt > 0.75 GeV/c/nucleon at which the yields start to drop. 
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Fig. 7. Elliptic flow parameters i>2 for neutrons (triangles) and hydrogen isotopes (stars, top panel) 
and their ratio «2,n/f2,h (bottom panel, dots) for central (b < 7.5 fm) collisions of 197 Au+ 197 Au 
at 400 MeV/nucleon as a function of the transverse momentum per nucleon pt/A in comparison 
with the UrQMD predictions for 7 = 1.5 (a-stiff) and 7 = 0.5 (a-soft) represented by the dashed 
lines. 

The results obtained for the second Fourier coefficient V2 describing elliptic flow 
are shown in Fig. [7] Here the combined data set for central and mid-peripheral 
collisions (PM3 to PM5) is used and compared to the corresponding UrQMD pre- 
dictions for b < 7.5 fm. The global rise of the absolute magnitude of v 2 with p t 
is well described even though the approximately 15% correction for the dispersion 
of the reaction plane (cf. Ref!^|) is not applied to the data. In contrast to the 
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directed flow and as expected from Fig. [5l the squeeze-out of neutrons depends 
significantly on the symmetry term chosen for the calculations. It is considerably 
weaker in the asy-soft case while the predictions for the asy-stiff case and for the 
hydrogen isotopes practically coincide (Fig. [71 upper panel). 

For a quantitative evaluation, the ratio of neutron-over-hydrogen flows is pro- 
posed as an observable which should be insensitive to uncertainties resulting from 
the experimental determination of the reaction plane and from the matching of the 
experimental and theoretical impact-parameter intervals. The comparison shows 
that the experimental ratios scatter in-between the predictions for the asy-soft and 
asy-stiff cases (Fig.[7j lower panel). A linear interpolation between these predictions, 
averaged over 0.3 < Pt/A < 1.0 GeV/c, yields 7 = 0.94±0.21. A smaller but within 
errors consistent value 7 = 0.52 ± 0.30 is obtained if the comparison is restricted to 
mid-peripheral impact-parameters 5.5 < b < 7.5 fmpS The power law coefficients 
depend weakly on the symmetry energy at saturation chosen in the parameteriza- 
tion. The slight increase of 7 if this value is lowered confirms that densities above 
saturation are probed. Other systematic uncertainties have been found to remain 
within A7 sa 0.2. Together with the kinetic term proportional to (p/po) 2 ^ 3 , the 
squeeze-out data thus indicate a moderately soft behavior of the symmetry energy 
at supra-saturation densities. 

This result can be considered as, within errors, consistent with the density de- 
pendence deduced from fragmentation experiments probing nuclear matter near or 
below saturation^ and with the slightly softer density dependence resulting from 
the analysis of the pygmy dipole resonance in heavy nuclei!^ It also shows that the 
super-soft density dependence resulting from the IBUU analysii^ of the tt~ /tt + 
yield ratios can, at this time, not be considered as a unique consequence of experi- 
ments probing higher densities. Pion yields are expected to be produced during the 
high-density phase of the reaction^! but are also subject to considerable in- medium 
effects which have to be controlled with high precision!^ More work, possibly in- 
cluding additional observables and focussing on the consistency among them, will 
be needed in order to arrive at firm conclusions regarding the symmetry energy 
at high density. In view of the far-reaching and mani-fold consequences, this task 
deserves highest priority. 

Illuminating discussions with Lie-Wen Chen, Bao-An Li, M. Di Toro, and 
H.H. Wolter are gratefully acknowledged. This work has been supported by the 
European Community under contract No. HPRI-CT-1999-00001 and FP7-227431 
(HadronPhysics2) . 
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